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We have recently introduced an approach for studying perturbatively classical and quantum 
canonical general relativity. The perturbative technique appears to preserve many of the attractive 
features of the non-perturbative quantization approach based on Ashtekar's new variables and spin 
networks. With this approach one can find perturbatively classical observables (quantities that have 
vanishing Poisson brackets with the constraints) and quantum states (states that are annihilated 
by the quantum constraints). The relative ease with which the technique appears to deal with 
these traditionally hard problems opens several questions about how relevant the results produced 
can possibly be. Among the questions is the issue of how useful are results for large values of the 
cosmological constant and how the approach can deal with several pathologies that are expected to be 
present in the canonical approach to quantum gravity. With the aim of clarifying these points, and to 
make our construction as explicit as possible, we study its application in several simple models. We 
consider Bianchi cosmologies, the asymmetric top, the coupled harmonic oscillators with constant 
energy density and a simple quantum mechanical system with two Hamiltonian constraints. We find 
that the technique satisfactorily deals with the pathologies of these models and offers promise for 
finding (at least some) results even for small values of the cosmological constant. Finally, we briefly 
sketch how the method would operate in the full four dimensional quantum general relativity case. 

I. INTRODUCTION 
A. The problem 

In the last few years, important developments in canonical quantum gravity have taken place, culminating with 
the recent formulation of two mathematically consistent , non-trivial canonical quantizations of general relativity 
in terms of Ashtekar new variables and spin networks. Instrumental in these developments have been the underlying 
advances in mathematical techniques for dealing with infinite dimensional nonlinear spaces, like the theory of cylin- 
drical functions and associated measures || , and the introduction of spin networks to eliminate the over-completeness 
of the Wilson loop basis Q (for a recent summary, see the review article by Rovelli || ) . In spite of the advances, 
most of the results obtained from these theories up to now concern statements made at a kinematical level (without 
imposing the Hamiltonian constraint). 

In order to start discussing if these theories contain the correct semiclassical physics and are therefore physically 
viable theories of quantum gravity, we need methods to introduce semiclassical states (encouraging recent kinematical 
results can be seen in |(|) but also to probe their dynamics. In particular a current open issue is if the semiclassical 
physics can be discussed at a kinematical level only or requires the imposition of the correct dynamics. Probing the 
dynamics has always been a problem in canonical general relativity . In any gauge theory, the only physically relevant 
quantities are those that are invariant under the gauge symmetries of the theory, or in the canonical language, that 
have vanishing Poisson brackets with the constraints that represent the gauge symmetries. At a quantum mechanical 
level, the physical states are those that are annihilated by the constraints. In general relativity we do not have any 
example of a quantity that classically has vanishing Poisson brackets with the constraints at least in the case of a 
compact manifold. In fact, it is strongly suspected that such quantities will never be constructed in closed form 
Finding quantum states that are annihilated by the constraints is also a challenge. 



x An exception to this may be the quantity constructed in fel, corresponding to the holonomy group of the spin connection. 
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B. An approximation scheme 



Most researchers in this area do not consider these problems a fundamental obstruction, since in order to do mean- 
ingful physics one only needs approximate expressions for the observables, as was advocated long ago by Bergmann 
and Newman JlC| ] . Finding an approximation method that does not destroy the non-perturbative nature of the canon- 
ical treatment is, however, a challenge. We have recently proposed one such method ItOj . It is based on considering 
general relativity coupled to a cosmological constant and taking the A — > oo limit p2[ . The construction goes as 
follows: 

For general relativity with a cosmological constant, the Hamiltonian constraint reads, 



H(N) = H A=Q {N) + A / d 3 xN(x) v / detq(x) (1) 



where N(x) is a smearing function, H\ = q(N) is the Hamiltonian constraint without a cosmological constant and detq 
is the determinant of the spatial metric. If one now considers the limit A — > oo, and re-scales the constraint by 1/A, 
one is left with a theory, which we will call "zeroth order" theory for which the Hamiltonian constraint is, 

H^{N) = J d 3 xN(x) v / detq(x), (2) 

that is, the Hamiltonian constraint is just the square root of the determinant of the spatial metric. In addition to this 
the theory has the ordinary diffcomorphism constraint and if one uses Ashtekar variables there will also be a Gauss 
law constraint (both are independent of A) . Imposing classically the Hamiltonian constraint of the zeroth order 
theory, one immediately sees that it corresponds to metrics of identically vanishing determinant, that is, degenerate 
metrics. 

It is easy to construct quantities that have vanishing Poisson bracket with the Hamiltonian constraint of the 
zeroth order theory. For instance, one can consider any function depending only on the three metric (or if one is 
using Ashtekar variables {Ef , A J b ), the densitized triads Ef) and not on its canonically conjugate momenta. To have a 
genuine observable the quantity should also have vanishing Poisson bracket with the diffeomorphism constraint. There 
are standard ways of achieving this coupling the theory to matter (l^J3 . A novel point is that in our approach, since 
the extra matter couplings in the Hamiltonian are higher order in A - these techniques yield genuine observables for 
the zeroth order theory, whereas normally they just construct quantities that have vanishing Poisson bracket only 
with the diffcormophism constraint but not with the Hamiltonian constraint. Since we will not need this technique 
for the examples we consider in this paper, we do not include further details here (see a discussion in p]]). 

If one assumes that the observables of the theory are power series in the inverse cosmological constant, 

Oa(it,E) =O [0 \n 1 E)+k- 1 O {1) {ir,E) + ..., (3) 

and one requests that these observables have vanishing Poisson brackets with the Hamiltonian of the theory (more 
precisely, one would like that the Poisson bracket be proportional to the constraints of the theory, for simplicity we 
will just demand that they vanish, but it is immediate to extend the construction to the case in which they vanish 
on-shell)and we expand such requirement in powers of A -1 one gets 

|o(0),ij(0)} =0 (4) 
[0^,H^)+{0^,H^) =0. (5) 

The first equation determines and the second one leads to a linear partial differential equation for The 
construction can be readily continued to higher orders. In all cases one obtains a linear partial differential equation, 
albeit with a more and more complex inhomogeneous term. It should be noted that one can obtain many observables 
starting with different solutions to the first equation. The linear partial differential equations are not hard to solve, 
given that the coefficients of the derivatives are functions of the triads, whereas the derivatives are with respect to the 
connections. Given the simplicity of it, the system is always integrable (although the solutions might be pathological, 
as we see in the examples later in the paper) and therefore yields all the observables of the theory. 

This amazing simplicity also has a quantum counterpart, at the time of finding states that are annihilated by the 
Hamiltonian constraint. If we now focus on the formulation of the quantum constraints in terms of the spin network 
representation of quantum gravity, for instance as discussed by Thiemann jjj, one starts by considering our zeroth 
order Hamiltonian, which is closely related to the well-understood volume operator. It is easy to find eigenstates for 



2 



this operator, in particular ones with vanishing eigenvalues. When one considers the first order corrections, because 
the zeroth and first order Hamiltonians have a well defined algebraic action in the space of spin network states, one 
gets a completely manageable (albeit complicated) system of equations. We will not go into details of the quantum 
theory of full general relativity in this paper, since it would make the paper much less readable as we would have to 
introduce quite a bit of (readily available but quite detailed) technology of spin networks. We will only present, for 
completeness, a brief of what the issues in the full quantum theory are in the last section of the paper. 

Let us however describe in general lines how the quantum perturbative approach would be for a generic system 
with a single Hamiltonian constraint (we will discuss the case of more than one Hamiltonian as a model of what can 
happen in a field theory in a subsequent section) . 

Consider a system where the quantum Hamiltonian constraint is of the form H (A) = + XH^ where A is a 
parameter that we consider small. We then consider quantum states that admit a power series expansion in powers 
of A, 

|V(A) >= |^ (0) > +A|V^ (1) > +X 2 \^ 2) >+... (6) 

What one wants to do now is to solve the theory perturbatively. This requires some comments. What does it 
mean to solve perturbatively a quantum mechanical constrained system? The approach we will take is the following: 
consider the eigenvalue equation for the Hamiltonian constraint, 

H(\M(\) >=E(\)\i>> (A), (7) 

and consider its solution up to a certain order in A, finding the zero-energy eigenstate within such a theory. That is, 
what we end up doing is bound-state perturbation theory (provided the spectrum of is discrete) and finding one 
set of particular states, those with zero eigenvalue. A good example to bear in mind about the method is to consider 
an atom in a (weak) magnetic field and think that one is looking for states with a given energy. If one goes to a 
certain order in perturbation theory in terms of the magnetic field, one will be able to achieve the given energy level 
only for certain values of the magnetic field. That is, in this approach the perturbative parameter is quantized. It is 
also clear that the perturbed states that achieve a certain energy level for a given value of the perturbative parameter 
do not necessarily stem from perturbing the zeroth order theory level of that particular energy. It is important to 
notice this in the gravitational case, where one looks for states with vanishing eigenvalue. The perturbed states of 
interest will in general not have H^\^ ^ >= 0. 

So we have cast the problem of solving quantum mechanically a constrained system as finding perturbatively an 
eigenstate for a given Hamiltonian in a quantum mechanical system. This problem is well defined and, as we will see 
in the examples, solvable. For the full general relativity case, the zeroth order Hamiltonian is quite close to the volume 
operator, so we do have a discrete spectrum to start with. However, because it is a field theory, one really has infinitely 
many Hamiltonians and certain consistency requirements have to be met. Therefore the quantum mechanical analogy 
we made may turn out to be too naive. The last example we consider in this paper is geared towards clarifying that 
point. 

Why should one care about a regime in which A is large in Planck units, given that the current observables value 
in such units is 10~ 132 ?. Our attitude is the following: we have now at least two candidates for theories of quantum 
gravity. We need to probe them in any regime we can get a handle on in order to help decide if the theories are 
plausible or not from a physical stand point. If we can do this, we can then see in detail what regime of A we can 
really access. Perhaps it is a bit better than what one initially expects and the results can be of interest to early 
universe cosmology. Or perhaps one is lucky enough to encounter a situation similar to other power series expansions 
(like the strong coupling limit in the lattice or the large N expansion in QCD) where the results are robust and useful 
way outside the naively expected regime. We will see some hints of this in one of the examples of this paper. 

C. Potential difficulties of the method and layout of answers provided in this paper 

The above observations make the whole approach worthwhile exploring further. However, the surprising way in 
which it appears to deal in a simple way with questions that in the full theory seem unapproachable raises a certain 
level of skepticism. Several pointed questions can be asked that make the whole approach look questionable. Among 
them, 

a) Does the approximation have any hope of working for values of A smaller than the Planck scale? 

b) General relativity is a complex enough dynamical system, that some pathologies might be expected. For instance, 
solutions might exhibit some level of chaos. It could also happen that other pathologies appear, like the full theory 
having less observables than the zeroth order one. Docs our method always produce 2n — 2m observables (where n 
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is the number of degrees of freedom and m the number of constraints)? With a method that reduces everything to 
linear ODE's, is there not a risk of "papering over" these subtleties. 

c) What sort of approximation is one getting at a quantum level? What is meant by an approximate state, or an 
approximate solution to a quantum constraint? 

d) When one applies the quantum perturbative approach in a field theory, one effectively is dealing with infinitely 
many Hamiltonian constraints. Is the method still viable or one gets an inconsistent set of conditions on the perturbed 
states, leading to an empty theory? 

e) The zeroth order Hamiltonian constraints we consider seem to admit many more quantum states than the full 
Hamiltonian constraints. Does the method produce spurious states for the full theory or correctly notices the further 
limitations that should appear at higher orders? 

This paper attempts to provide (at least partial) answers to these questions by probing the perturbative method in 
model systems. In section II we will consider the application of the technique to Bianchi models. This will allow us to 
probe questions a) and b). From Bianchi models we will learn that, surprisingly, when one evaluates the approximate 
power series that we get for the observables on solutions to the equations of motion, the cosmological constant drops 
off from the expressions. Evaluated on the solutions to the equations of motion, the expression for the observables 
become a sum of terms that approximate better and better a constant of motion without reference to the cosmological 
constant. For a Bianchi I recollapsing universe we will explicitly demonstrate that with only a few terms one finds an 
expression that maintains an almost constant (less than 10% variation) value across 80% of the life of the universe. 
We will also see that in the case of Bianchi IX one can easily construct observables, but it is likely that they will 
remain good approximations only for short periods of time when one enters the chaotic domain. We will also see that 
one does not obtain more observables than the required ones. One can obtain many apparently different power series 
with the method but they just represent reparametrizations of the basic observables. We will also learn how to deal 
with possible singularities in the expressions obtained for the observables as functions of the phase space. 

In section III we will discuss the heavy asymmetric top in three dimensions. This is a Hamiltonian system that can 
be viewed as a pathological constrained system if one considers it for a fixed value of the energy. It is pathological in 
that it has less observables than four. We will see that our method only apparently generates too many observables. 
On close inspection we will find that several of the expressions are really not good candidates for observables because 
they are expressions that have support only on limited regions of phase space, that are in general not preserved by 
evolution. In the end, the method produces the correct number of observables for the model. 

In section IV we explore a model consisting of two harmonic oscillators with constant energy difference. This system 
has received quite a bit of attention, for instance see Jl5| . Unless the ratio of the frequencies of the two oscillators 
is a rational number, the system is pathological in that it only admits one quantum state. We will show that our 
perturbative approach indeed notices this fact and produces correct approximations for the quantum states in all 
cases. 

In section V we explore a model of coupled harmonic oscillators with two "Hamiltonian" constraints. General 
relativity being a field theory one has infinitely many Hamiltonian constraints (one per spatial point). Therefore 
the other models we consider in the paper seem a bit short in mimicking this feature, since they are mechanical 
systems with only one constraint. In this section we increase modestly this feature by considering a model with two 
Hamiltonians. We will see that the method deals with no problem with such a system. 

We will end with a short reflection on the implications of the results of this paper and what challenges remain for 
the application of the approximate scheme in the context of full classical and quantum general relativity. 

II. BIANCHI MODELS 

As a first example of the application of the technique, we will discuss the issue of finding observables for Bianchi 
models. Bianchi models are homogeneous cosmologies, and as a consequence the Einstein equations become ordinary 
differential equations. This allows to achieve, in some cases, quite a bit of progress towards their solution. 

Following (l6|| , we represent the Bianchi models using as variables the (diagonal) components of the triads E l and 
the canonically conjugate momenta Aj. All class- A Bianchi models admit a simple Hamiltonian formulation in terms 
of these variables, but we will concentrate our attention in two models, the Bianchi I model and the Bianchi IX 
model. The former is the simplest model and the latter the one with the richest dynamics. Bianchi I models start 
from an initial singularity and expand anisotropically. Depending on the sign of the cosmological constant they may 
expand forever or recollapse. Bianchi IX models approach the singularity through an infinite sequence of oscillations 
of increasing complexity. 

In terms of the Ashtekar variables the Bianchi I and IX models with a cosmological constant are described by a 
Hamiltonian constraint, 
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H = E 1 E 2 E 3 - i (E 1 E 2 (A 1 A 2 - eA 3 ) + E 1 E*{A X A 3 - eA 2 ) + E 2 E 3 {A 2 A 3 - eA x )) (8) 

where e = 1 corresponds to the Bianchi IX model and e = to the Bianchi I model. The Gauss law and diffeomorphism 
constraints are identically satisfied. 

To apply our construction we need to start by choosing an observable for the zeroth order theory. In general such 
quantities are given by = F(E 1 ,E 2 ,E 3 , E 1 A\ — E 2 A 2 ,E 1 A X — E 3 A 3 ). If one sets out to find four independent 
observables one could start, for instance, with the following independent choices for zeroth order observables, 

(9) 
(10) 

(11) 
(12) 

Let us now concentrate on the study of the Bianchi I cosmology. In this case O34 are both exact observables for the 
full model with a cosmological constant. Starting with of one finds the following first order correction, 

0« = (-E^E 2 A 2 + (E 1 ) 2 A 1 E 1 E 3 A 3 ) , (13) 

and in addition to this one has the solution to the homogeneous equation, which is given by an arbitrary function 
that commutes with H^, that is, similar to 0^°\ The expression for the Bianchi IX model is very similar to the one 
we just introduced, it has an extra term —eA x . 

In order to investigate the behavior of the approximate observable we are considering we will evaluate it on an 
exact solution of the equations of motion, the metric of reference jlTj , 

ds 2 = -dt 2 + V{t) 2 ' 3 \?,(t)i cos ^dx 2 + Z(t)? co <~<-^)dy 2 + £(t)s cos (T+Hdz 2 l (14) 



of 


= E 1 




of 


= E 2 




of 


= E 1 A 1 


~E 2 A 2 


of 


= E 1 A X 


-E 3 A 3 



where, 



V(t) = ^ (15) 



= 2(1 -cos art) 
uj sin Lot 

and where lo = \J — 3A and where A < 0. This solution corresponds to a Bianchi I universe that expands out of a Big 
Bang at t = but after a while the cosmological constant leads it to recollapse at a time Lotp — ir. The Ashtekar 
variables for this metric (particularized to 7 = 0) read, 

E 1 = (cos(l/2u> t)f 3 (17) 

e2 = E 3 = ff2^tan(l/2a;f) (sin(o; t)f 3 

to 

#2(2 sin(u; t) (cos(l/2 cot)) — 3 sin(w t) — 4 sin(l/2 Lot) cos(u; t) cos(l/2 lo t) ) 

M = 1/6 ^ jj, - 2 '- (19) 

^/sm(ut) (sin(l / 2 Lot)) 2/ 3 (cos(l/2 Lot)) 2 

A 2 =A 3 = l/3 MV 2 "*)". (20) 
' ?Jcos{l/2Lot) V ^ 

As we can see, for the solution in question E 2 and E 3 vanish at t = 0. If we look at the expression for the first 
order correction we found, this implies that it diverges for t = 0. This is an example of what we stated earlier, namely 
that the expressions for the observables we found can sometimes be singular in certain points of phase space. One 
can construct a first order correction that is well behaved at t = by making use of the free function F^ that solves 
the homogeneous equation. If one adds to the first order correction the solution to the homogeneous equation given 
by (E 1 A 1 - E 2 A 2 )(E 1 A 1 - E 3 A 3 )(E 2 E 3 )~ 1 , it becomes, 

(1) = A 2 A 3 (21) 
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and this quantity is obviously well behaved at t = 0. 

If one now proceeds to compute the second order correction, one finds a similar situation with 
respect to the divergence at i = 0. This can be corrected by adding the function F 2 = 
-1/3 {{E l A x - E 2 A 2 ){E 1 A l - E 3 A 3 )) 2 (E 2 E 3 )- 2 (E 1 )- 1 and one obtains for the correction, 



(2) 



1 



? 3\2 



-2{E v y A{E Z E 6 A 2 A 3 + {E'f A{Ai,{E 3 ) 



6(E 2 E 3 ) 2 E 1 
-2E 1 A 1 E 2 {E 3 ) 2 A 2 {A 3 ) 2 + {E 1 ) 2 A 2 A 2 {E 2 ) 2 
-2E 1 A 1 (E 2 ) 2 E 3 A 2 2 A 3 - 2A\(E 2 ) 2 A\e\) 



(22) 
(23) 



and a similar subtraction yields the third order correction 
1 



(3) = _ 



{-%{E x f A 3 (E 2 ) 2 : E 3 A 2 A 3 - GiE^A^iE^A^ 2 - 12E 1 A^E 2 ) 3 \E 3 ) 2 A 3 A 2 (24) 



SO^ 2 ^ 3 ) 3 ^ 1 ) 2 

-UE^^fiE^AlA 3 - ^(E 1 ) 2 A\(E 2 ) 2 (E 3 ) 2 A 2 A\ - 30E 1 ) 2 A 2 E 2 (E 3 ) 3 A 2 A 3 
-3(E 1 ) 2 A 2 (E 2 ) 3 E 3 A 3 A 3 ~ ^(E 1 ) 4 AfA 3 E 3 A 2 E 2 + 2(E 1 f A\A 2 {E 3 ) 2 + 2{E 1 fA\(E 2 ) 2 A\ 
-6A 3 (E 2 ) 3 A 3 (E 3 ) 3 + SAKE^iE^A 3 + G^AfA 3 ^ 3 ) 3 ) 

If we now evaluate the above expressions for the exact solution we discussed above one gets (x = ujt) , 



01' 



- (0) = cos 4 / 3 



(1) 



A 



= 1/3 



2cos 2 (§) + l 
cos 2 / 3 (f) 



01 



(2) 



A 2 



1/27 



-19cos 4 (|) + 19 cos 4 (|) cos 2 (x) - 6cos 2 (|) + 6 cos 2 (|) cos 2 (x) 



—4 cos( — ) sin(x) cos(ir) sin(— ) — 2 + 4cos 3 ( — ) sin(x) cos(x) sin( — ) + 2 cos 2 (x) cos 8 / 3 ( — ) (— 1 + cos 2 (x)) 

Z Z Z Z -I L £ 



(25) 
(26) 

(27) 
1 



(3) 



A 3 



1 

810 



-568cos 6 (^) - 15+ 15cos 2 (a;) - 194 cos 4 (^) + 488 cos 6 (^) cos 2 (x) + 354cos 4 (^) cos 2 (x) 



— 47cos 2 ( — ) cos 2 (x) — 33 cos 2 (— ) — 16cos 3 (— ) sin(x) cos(x) sin( — ) — 72 cos(— ) sin(x) cos(x) sin(— ) 

Z Z Z Z Z z 

+88cos 5 (|)sin(x)cos(x)sin(|) cos 14/3 (|) (-l + cos 2 (x)) 



(28) 



The first observation is that the above expressions for the correction have lost the dependence on A stemming from 
the perturbative approach (there is a trivially re-scalable dependence on A through x = u)t). Therefore, if the above 
expressions converge, the convergence will be independent of A. This can be seen in the following plot, in which we 
show the observables as functions of x. We see that in spite of the lack of dependence in A, the curves do converge, 
the third order one maintaining an almost constant value for the longest period of time (more than half the lifetime 
of the universe). The various corrections are not well behaved near the big crunch. 



6 



0.5 



Zeroth order 


/ 1 
/ / 


First order 


/ / 


Second order 


/ /' 


Third order 


1 
1 




i 

i 

o / 1 


\ * — 


2 V 






s 

s 


— ~ ■ - \ 


\ 


1 \ 


\ 







\ \ 
\ \ 




\ 




\ 




\ 




\ 




\ I 



FIG. 1. The perturbative observable evaluated for a Bianchi I cosmology with negative cosmological constant. The time 
runs from zero, the Big Bang, to tv, the Big Crunch. The four curves correspond to better and better approximations. 



It should be remembered that the above calculations were particularized for 7 = 0. The value of the observable will 
in general be 7-dependent. One way of understanding this is that the "physical" information content of the observable 
is to give a measure of the 7 parameter, which is associated to the anisotropy of the model. For instance, the exact 
observable O3 = A 1 E 1 — A 2 E 2 takes the value O3 — — cos 7 + sin 7. Since the solution to the equations of motion 
we are considering is parameterized by a single constant, all observables become given functions of that constant. 

As was shown in the above calculations, the Bianchi IX model appears as only slightly more complex, and similarly 
any class- A Bianchi model can be treated in the same way. Because of the local nature of our calculations (in time) , 
it is not surprising to find observables in the Bianchi IX case. In that cosmology, because of the chaotic oscillations 
the dynamics undergoes, it is unlikely one will be able to find expressions for the observables that are global, but 
expressions like the ones we find should face no difficulty. 

As we saw, the zeroth order observable one starts from has a huge degree of ambiguity. There is also a lot of 
ambiguity in each perturbative step, due to the presence of the solution of the homogeneous equation. It is evident 
that certain choices will lead to more useful (or in some cases shorter or even polynomial) expressions, as we exhibited 
when we manipulated the homogeneous solutions to avoid expressions that are singular at certain points in phase 
space. 

Summarizing, we learn the following lessons from the Bianchi example we have just discussed. First of all, we see 
that the approximation for the observables we get, when evaluated on a solution of the equations of motion, becomes 
independent of the cosmological constant. This suggests that the perturbative approach will remain good even when 
the constant is small. In fact, we see in the example, that the observable is well approximated for a significant fraction 
of the lifetime of the universe in question. We also see that the presence of chaos in the Bianchi IX model does not 
preclude us from finding approximate expressions for the observable, though it is likely that it will impose further 
restrictions on the domain of validity (in terms of evolution time) of the approximate observable. It is possible to test 
these hypotheses with numerical evolutions of the Bianchi IX cosmology. 



III. THE ASYMMETRIC HEAVY TOP WITH A FIXED ENERGY VALUE 

As another example of the application of the perturbative Hamiltonian ideas we are proposing, we analyze a model 
which in spite of its simplicity, adds pathological behaviors that are interesting to analyze under the proposed scheme. 
The model is an asymmetric heavy top with constant energy. This system only has one observable, the momentum 
canonically conjugate to the angle around the vertical axis. We will assume that the asymmetry in the top is small 
and treat this model as a perturbation of a symmetric heavy top, which has two observables. Will our technique 
notice that one loses an observable when one loses the symmetry? 

Consider a heavy symmetric top with moments of inertia I\ — I2 and 1%. Now suppose one sticks on two small 
masses at opposite sides of the surface of the top, forming a line perpendicular to and intersecting the top's third 
axis at its center of mass. Then the moments of inertia along directions 1 and 2 will be slightly different. Let us 
quantify this as I 2 = h +A, with A a small perturbation (we use the lowercase to avoid confusion with the cosmological 
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constant, although the parameter plays the same perturbative role in this problem as the inverse cosmological constant 
in full general relativity). 

The Lagrangian for such a system will be given, in terms of the traditional Euler angles (see for instance [Q), 

L = ]-h(e 2 + tp 2 sin 2 e) + ^ I 3 (ip 2 + ip 2 cos e) 2 + ^(tp sin 9 cos V - sin ip) ^ - cos 6 (29) 



2 V r / 2 V r / 2 

where £ is the length along the top third axis from the top base up to the center of mass. To introduce a canonical 
formulation, we compute the canonical momenta, 

dL 

pg = — r = 1 16 + \9 sin 2 ip — \ip sin 9 cos if) sin ^ (30) 
89 

and similar equations for the other momenta, which yields an invertible linear system of equations from where one 
can compute the velocities as functions of the canonical momenta. One can then perform a Legendre transform and 
obtain a Hamiltonian H{pg,p Vl p^, 9, ip, A). For small values of the perturbative parameter, one can expand, 

OH 

H = H(\ = 0) + \— (X = 0) = H + XH U (31) 
uX 



where the general form of the Hamiltonians is, 

ffo^+E^ftPj, (32) 



2Ji 

3 

H 1 = =) j d i:j (9, if)p l p j (33) 

with pi — pe,P2 — Pip,P3 — Pip- We would like to treat the system as an example of constrained system. Therefore 
we request that the energy have a given value H = Eq, that is, 

H (o) _ E (o) + XH (D = Q (34) 

Before treating the case of the top with constant energy, let us first recall how one formulates the usual asymmetric 
top, but as a parameterized system. This allows a more natural contact with the constrained systems we will handle 
later. The parameterized top is obtained by introduce a time q° and its canonically conjugate momentum p$, such 
that, 

Po-H( Pu q l ) -0 (35) 
where q 1 = 9,q 2 = ip, q 3 = ip. The observables are defined by, 

{O lP o-H( Pl ,q*)} =0, (36) 

which yields, 

80__dOdI£ + dOdH__ Q 
dq° dq l dpi dpi dq % 

This system is solved by finding the solution to the following set of differential equations, 

ctf =i dtf = _dH_ dp^ = dH_ 
ds ' ds dpi ' ds dq l 

The first of these equations can be integrated immediately, q° = s + C. The other two equations are just Hamilton's 
equations of motion, which integrated yield, 

<f = f i (Pi,4,s),Pi=9i(Pi,qls), (39) 
if one substitutes s by — q° in the above expressions, one obtains the six constants of motion, 
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3? = /i(Pi,? i ,-g°) > pi = 9i(Pi,q i ,-q°), (40) 

which are observables of the parameterized system. A crucial element in the solution of this system, that allows to 
obtain the six constants of motion explicitly is that the solution of the equation for s in terms of q° is straightforward. 
We will see that this poses problems in the next example. 

Let us now study the case of interest, the (slightly) asymmetric heavy top, with Hamiltonian constraint ( |34"| ) and 
assuming observables of the form O = Oq + XO\ + . . .. Let us start with the equations to zeroth order in A, 



{o (0) , 



ff(o)_£(o)j =0j (41 ) 



with 



H (o ) = (^_ + j2 Cij{e)PiPjj (42) 



which leads to the PDE, 



pi dO A dO Q dO ddj 



which again can be solved by the substitution technique, 

*,J'=2 

At zeroth order the system has two trivial constants of motion, 

p i =p° i , i = 2,3. (46) 

One can in principle integrate the remaining equations to yield, 

= 0(Pi,qis),Pe=Pe(plqls) (47) 
<p = <p(Pi,Qo,s)^ = il>(pi,Q(„s), i = l,3. (48) 

As before, one can (formally) invert the solutions to obtain the constants of motion as functions of phase space by 
solving for s and substituting back in the following equations: 

<7o = MPi,<l\s),Pi = 9i(Pi,<l\s), i = l,3 (49) 
p i =p° i , i = 2,3. (50) 

We therefore appear to have found again six constants of motion as in the case we studied before. However, the above 
equations involve periodic functions with arguments with incommensurate periods. One will only determine s up to 
multiples of the period and when substituting in the other equations will end up with multivalued functions. If one 
chooses a particular branch, the resulting observables only take values on a portion of the accessible phase space. 
They are therefore not acceptable as observables for the system. 

This is an important observation, since we will see that this is the mechanism which prevents our technique from 
generating observables for systems that do not have them. 

Let us now consider the first order corrections. The equation to be solved are, 

{O 1 ,H -E } + {O ,H 1 }=0 (51) 
and if we start from Oq = p^p, we are left with, 

Pi dOi A dO x dOi ^ , , -A ddij 

1 ^ i,j=2 ^ FV i,j=2 i,j=2 Y 
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Now the PDE is non-homogeneous. We can still solve it with our parameterization technique via, 



= H * = 2,3 (54) 

= XI 4*wj> * = 2 ' 3 ( 55 ) 

$1 = (56) 
as 

& = E— . = 1,3. (57) 

As before, we will need to solve for s and replace it in the equations to obtain expressions for the initial values as 
functions of phase-space. We do run, however, into the same difficulty as before, we cannot obtain a uniquely defined 
expression for s and therefore we cannot construct a first order correction for the observable chosen. The procedure 
therefore fails to work. 

What happens if we choose Oq — p v instead? In that case {Oq, Hi} — and therefore the equation for 0\ indeed 
has a solution. We will therefore recover, order by order, an expression for the single observable the system admits. 

The general conclusion seems to be that, although the linear PDE's of our procedure in principle can always be 
solved, requiring that the solution be a well defined function on phase space (and not only on a portion of it) severely 
restricts the observables that can be obtained by the technique, and apparently prevents the technique from generating 
spurious observables in systems with less observables than degrees of freedom. 



IV. TWO OSCILLATORS WITH CONSTANT ENERGY DIFFERENCE 



A. The model 



Let us now consider examples of the application of the method at a quantum level. Haji'cek and others ]15[] have 
considered a model of constrained system, consisting of two harmonic oscillators with a fixed energy difference. The 
possible eigenstates of the system therefore consist (if one uses the standard inner product for the oscillators) in 
situations in which oscillators have the appropriate quanta of energy to satisfy the constraint that the difference of 
their energies be constant. It is clear that this can only happen if the ratio of the frequencies of the oscillators is 
rational. Otherwise, the system would admit at most one quantum state. This is the pathology of this system: if 
the ratio of the frequencies is not rational, the quantum system only has one state satisfying the energy constraint. 
If it is rational, the system admits infinitely many states. One could therefore consider starting with a zeroth order 
model consisting of two coupled oscillators with rational ratio of frequencies and perturb the system by adding the 
Hamiltonian of two oscillators with an irrational ratio. The zeroth order system has infinitely many quantum states 
and the perturbed system has only one state. Can our method detect and correctly handle this situation? This is 
what we attempt to probe in this section. 



B. Zeroth order Hamiltonian 



The Hamiltonian constraint for the system in question simply reads, 

H = \{p\ + wlxl) -\{pl+ ulxi) + 1 = (58) 

where we have chosen the energy difference equal to one. The quantization of this model, viewed as a constrained 
system, has order ambiguities. Consider the transformation 

IT^i, Pi = \ —Pi, (59) 
n V 
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the constraint reads, 



H = h*J\{& + x\) - + xj) + 1 = (60) 

and introducing the complex variables, 

a 3 = 7^(% + a j = "^(% - ( 61 ) 

the constraint takes the form, 

= Tujj\ot\a.\ — Tuj}ia*2tt2 + 1 = 0. (62) 
We can then proceed to the quantization, 

Xj^Xj, pj^Pj (63) 

[xj,P k ] = i$jk (64) 

ctj — + &j, a* — ► aj, (65) 

and the quantum constraint can be written as iJ = _ffi — + 1 = 0, with, 

#j = ficjj (^fija^aj + (1 — Hj)cijo)^j (66) 

with ^ representing the ordering ambiguities. To obtain the quantum states that arc annihilated by this constraint, 
let us consider separately the spectra of both -ffj's. Since they are harmonic oscillators, we have, 

Hj = TiLUjUj + HKj, eigenvector \rij > (67) 

with K\ 2 arbitrary ordering-dependent constants. In order to solve the constraint, the values of n\,ri2 must satisfy, 

hio\ri\ + TiK\ — huj2n2 — hK% + 1 = 0. (68) 

To produce solutions, one can choose the constants Kj. For instance, for n\ = ri2 = 0, we get U(Ki — K2) + 1 = 0. 
In order to have more than one solution, one needs to have, 

bJ\Tl\ — U12U2 — 0, (69) 

or, 

^ = ^ (70) 

that is, the ratio of the frequencies of the oscillators must be a rational number. If 21 = £L then ri\ = Nr, n% = Mr 
with r an integer and the wavefunctions of the physical space are \Nr, Mr >= \Nr > \Mr >. 

C. Perturbed model: exact treatment 

We will now construct a model that is a perturbation of the above considered one but is still exactly solvable. We 
can then solve the model exactly and perturbatively and compare the results. Let us consider the system f| 

ui 2 

H = Hx-H 2 - x y x z + 1 = °- ( 71 ) 



2 Notice that this is equivalent, up to an overall rescaling, to add a perturbation corresponding to the energy difference of two 
oscillators, as we proposed doing in the introduction. 
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This system is equivalent to the original Hamiltonian constraint if we define ui' 2 = ^(1 + A) 1 / 2 . We can therefore 
readily solve it. For system to have more than one state, we have, 

LOinx + uj' 2 n 2 = 0, (72) 

or, 



ni _uj' 2 _ w 2 \/l + A 



n 2 ui\ u)\ 

ivinff fViof , 



therefore n\ = Pr and n 2 = Qr with r integer, and remembering that ^ = j^, in order to have a solution the 
perturbative parameter has to have a given value, 



P 2 M 2 

A = W~ L < 74 » 

In order to compare with the perturbative calculation we will perform in the next paragraphs, let us consider the 
condition for having more than one solution for this model, 



LOini — uj' 2 n 2 = = LUini — vl + Xuj 2 n 2 (75) 
which in the case A <C 1 reads, to first order, 

W in 1 -(l + ^) W2 n 2 = 0, (76) 



or, 

\ = 9 I 

.AT n 2 

In order to have Ac 1, the ratio n\jn 2 must be a rational number close to N/M. We can therefore write it as, 

m _ iVr + £ 
n 2 Mr + m ' 

with large r. Expanding, we therefore get, 



A = 2^I-1). (77) 



(78) 



A = 2 



2---2— + 0(r- 2 ). (79) 
iVr Mr y ' y ' 



Summarizing, the perturbed model can have infinitely many or a single quantum state, depending on the value of 
A. Let us now sec if treating the model perturbatively we recover the same result. 

D. Perturbative treatment of the model 

We consider the zeroth and first order Hamiltonians as, 

=Hi-H 2 + 1 (80) 
= -\ulx\ = (81) 

We start by considering the eigenvalue problem up to first order in A, 

(< + A0«|) (m + Aff«) = (e<°> + Ae«) (< O + A^|) . (82) 
Let us start by considering the resulting zeroth order equation, 

<^°>|ir< >=e< > <0 (o) |. (83) 
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The spectrum of is, 



e£°) = ?i ( Wl ni - u 2 n 2 + K 1 -K 2 ) + 1. (84) 



Here we need to distinguish if u>i/u>2 is rational or not. If it is rational, then the spectrum is degenerate and we need 
to apply ordinary degenerate bound state perturbation theory as described in any quantum mechanics textbook. If 
it is not rational, the spectrum is not degenerate. As we discussed in the introduction, let us concentrate on the case 
in which the frequencies are related by lo 2 /uji — N/M. The degeneracy of the spectrum can be directly seen in that 
it takes the same values for m = p + Nr, n 2 — q + Mr with given p, q and arbitrary values of r, that is, 

4°) = h (wip - to 2 q + Ki-K 2 ) + 1 (85) 

The eigenstates of H^ ' are, 

< <?WI =<p + Nr,q + Mr\ (86) 

in the number representation for the harmonic oscillators. We therefore find that the solution to zeroth order is given 
by (at the moment) an arbitrary combination of the given eigenstates, 

r 

Let us now consider the equation to first order in A, 

< < 4&\hM = 4°) < + e« < $)|. (88) 

Since (< <j>o\ + A < <f>i\) (\<po > +\\4>i >) = 1, to first order in A and we can choose the phase of |0i > to have, 

<0 O |0 O >=1, < ^ol^i >=< 01 1 0o >= 0Vp,o. (89) 

We now take the equation fl8Sj ) and project it onto the the eigenspace with eigenvalue e^) q . While doing this, the 
contribution from < cf>^l\H^ acting on this space cancels with the term e^l < (j>p~)q\ acting on the same space. We 
are therefore left with, 

< (f>^]j\H (1) \p + Nr,q + Mr >= e (1) < <pW\p + Nr, q + Mr >, (90) 
and recall that = l/2fajj 2 x 2 and in turn, 

(a 2 + of] a\ + (aX) + a 2 a\ + ala 2 a\ + (at) + 2a\a 2 + 1 

x\ = ± '— = ±-L = ^-L . (91) 

2 2 2 2 v ' 

Since the matrix elements of are between states of energy e^Jj, a 2 and (a^,) 2 do not contribute and the 

Hamiltonian is diagonal. Therefore of all the terms in the superposition defining < </>(°) | we are left with < p + Nr, q + 
Mr I, which normalizing yields, 

< 0o I =<p + Nr,q + Mr\ (92) 

and energy, 

= < 0o| 2 ^ a 2 2 + 1 |0o >= -\h^2 (q + Mr + . (93) 

Therefore the first order corrected energy is, 

e (°) + Ae« = hu 2 (^Lp - q - ±{q + Mr + + h {K x - K 2 ) + 1. (94) 

and as before we need to choose K\, K 2 in such a way that |0, > is a solution of the theory up to the order we are 
considering, i.e., K\ — K 2 + 1 — jHuj 2 = 0. If we want other solutions, we need to make the energy eigenvalue vanish, 
up to the order we are considering this implies, 
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A = 2 




(95) 



For small values of A we should recover the results of the exact calculation. We notice that A decreases with 
increasing r. In that limit we therefore have that, 



and this agrees with the expansion in small A of our exact calculation. 

It should be noted that the above calculation shows that starting from a given state in the zeroth order theory and 
choosing a perturbative parameter, we get a first order solution to the constraint. It is immediate to see however, 
that this can be accomplished by starting from an infinite number of sets in the zeroth order theory. Simply consider 
a state obtained by multiplying p, q and r times an integer. The same value of lambda in equation (^6|) ensures that 
the first order state is a solution. 

The calculation can be completed by considering the projection of (|8^) on the subspace of states orthogonal to the 
states of energy e£} q . Such projection determines the correction to the state, < <p^\. This can be straightforwardly 



The calculation can be repeated for irrational quotients of frequencies. In that case one is doing ordinary (non- 
degenerate) bound state perturbation theory. The calculations resemble very much the ones we did above, so we will 
not detail them here. Some comments are nevertheless in order. 

The zeroth order theory has only one state when the frequency ratio is irrational and many states when it is rational. 
In the perturbed theory, in both cases we can find values of A such that the theory has many states. That is, the 
perturbed model we constructed is such that depending on the value of the perturbative parameter one has either one 
or many states. If we do not choose the correct value of lambda (as determined by equation (Pq)), then the model 
only admits the state |0, > as solution. 

In particular, we have correctly addressed the situation posed in the introduction to this section: if the zeroth order 
model has rational frequencies (and therefore infinite solutions) but the peturbed model corresponds to a value of A 
that only admits one solution, we will not satisfy the energy constraint to first order and we are therefore left with the 
vacuum as the only solution (in the factor ordering chosen). Therefore the method seems to handle well the drastic 
reduction in number of states implied by the perturbation. 

It is suggestive that the model considered has a quantized value of A in the full treatment if one wishes to have 
infinitely many states. That behavior is reproduced correctly in the perturbative approach. However, the quantization 
of the perturbative parameter is a feature of the perturbative method as long as the zeroth order Hamiltonian has a 
discrete spectrum (which is the only case in which we can apply the method). It is interesting to notice that in the 
other models considered by Hajicek, the spectrum of the Hamiltonian is continuous and therefore we cannot apply 
our method, and the full models do not have any particular requirements on the frequencies of the oscillators. 

Another lesson from this example is that the lack of states in the perturbative treatment strongly suggests that the 
full model also lacks quantum states, unless there is some sort of non-analytical behavior in terms of the perturbative 
parameter. 



In the case of quantum gravity, one is dealing with a field theory, therefore one has an infinite number of constraints. 
In the spin network representation, if one considers networks with a finite number of vertices, the number of equations 
resulting from the constraints might be finite, but in realistic situations of semi-classical interest it will still be large. 
Therefore the quantum models we have considered up to now, which are quantum mechanical systems with one 
constraint, may not capture some of the aspects present in the gravitational case. It might be that dealing with a 
larger than one number of constraints implies further relations between the perturbative coefficients that make the 
system incompatible. It is somewhat unlikely that this will happen, since after all, assuming the full theory is solvable, 
if the solution admits a power series expansion in A, one should find a solution that is acceptable perturbatively to 
any order in perturbation theory. Nevertheless, it might be useful to consider a simplified model with more than one 
constraint to make sure the technique works. It is in particular interesting to see how, if one has several constraints 
to satisfy, one can achieve it with the same value of the perturbative parameter for all the equations. This is what 
we attempt in this section. 




(96) 



done. 



V. A MODEL WITH TWO HAMILTONIANS 
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A. The model 



Consider three different harmonic oscillators, 

Hi = ±(p? (97) 

and define the two "Hamiltonian constraints" , 

H^(1) = H 1 +H 2 -H 3 = (98) 
(2) = H 1 -H 2 + 2H 3 = 0. (99) 

The two constraints have vanishing Poisson brackets among themselves. We will have to request that the perturbed 
Hamiltonians also have vanishing Poisson brackets for all values of A and this will impose restrictions on our pertur- 
bative approach. 

As we discussed in the previous section, the eigenstates of the quantum version of the constraints are simply given 
in the number representation by |ni,n 2 ,n 3 >, 

i7 (0) (l)|ni,n 2 ,n 3 > = fi (wini + K x + uj 2 n 2 + K 2 - w 3 n 3 - K 3 ) \m, n 2 , n 3 > (100) 
H (0) (2)\n 1 ,n 2 ,n 3 > = h (wmi + K x - uJ 2 n 2 - K 2 + 2uj 3 n 3 + 2K 3 ) \n u n 2 ,n 3 > . (101) 

It is possible to choose the Ki in such a way that |0, 0, > has zero eigenvalue. The system admits other solutions 
only if 

MP , N 

^2 = ^^l' ^3 = (102) 

which leads to a system of equations for the n i; 

M P 

ni + —n 2 - -n 3 = (103) 

MP , s 

ni - ^n 2 + 2-n 3 = 0, (104) 

that has as general solution, 

ni = £PM, n 2 = -UPN, n 3 = -~2£QM, (105) 
with t an arbitrary integer. Therefore the physical states of the theory (states annihilated by the constraints) are, 

\£PM,-3£NP,-2£QM >, integer! (106) 

We will now consider a perturbative term that still allows the system to be solved exactly, 

H' 3 = H 3 + ±<44, (107) 

so, 

H(l) = H°(l)-^lxl (108) 

H(2) = H"(2) + ^lxl (109) 



The modification introduced by the perturbation is equivalent to changing lo 3 by u 3 = y/l + Xlu 3 . If we now consider 
the construction of quantum states \n\,n 2 ,n 3 > that are annihilated by both Hamiltonian constraints, and seek for 
values of rij that differ only slightly from the unperturbed ones, which is achieved by choosing p,q,s <§; t, 

m = £PM +p, n 2 = -WP£ + q, n 3 = -2£QM + s (110) 

we find that, 
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yr p p 

P+^q- (VTTX - l)-(-2£MQ) - VTTA-s = 



M 



P 



-P 



P - — q + 2(VTTA - \)-{-2lMQ) + 2VTT\-s = 0. 



(Ill) 
(112) 



The joint solution of these equations requires that p = Mb, q = —3Nb, and assuming that A is small, we find that, 

2Mb-^rS+^{2tMP-^rs) = Q (113) 

which implies that 

o Oh 1 

(114) 



s 2b 



IMQ IP y £ 2 ' 



B. Perturbative treatment 



We recall the form of the unperturbed Hamiltonian and the first order perturbations, 

HM(1) = H 1+ H 2 -H 3 , H^=- l -utxl 



H i -°\2) = H 1 +H 2 + 2H 3 , H^ = 2\ 



2 2 



with u)2 = %vi and uj 3 = ^Wi, and we wish to solve for quantum states, 

< | + A < | (h<® (a) + \HW (a)) = ( e (°) + Ae« ) (< <j><® | + A < 0« |) , 
The spectrum of H^°> is given by, 



= huji 



M P 
711 + N n2 Q n \ 

M P 
ni N n2 + 2 Q" 3 



+ h(K 1 + K 2 -K 3 ), 
+ h{K 1 -K 2 + 2K 3 ) . 



Once more, we have a degenerate spectrum, with, 

ni = iPM + p, n 2 = -MNP + q, n 3 = -21QM + s, 

so the energies are, 

£i 0) (P. 9. s ) = 



M P 



and the eigenstates are, 



< 1>?l J =< £PM + P. + g, -2£QM + s| 



M P ' 
P ~N q + 2 Q S , 



lh(K x +K 2 - K 3 ) 
+ h(K 1 -K 2 + 2K 3 ) , 



,p,q,s 



a = 1,2. 



(115) 

(116) 
(117) 



(118) 

(119) 
(120) 

(121) 

(122) 
(123) 

(124) 



Going to next order in perturbation theory, we again have to require (as in the example we discussed in the previous 
section) that the first order correction to the state satisfy, < <^°) l^ 1 ) >= 0. We now compute the first order correction 
to the eigenvalues, 



< 0(°> | ffW (a)\p, q, s, i >= e« < \p,q,s,£>, a = 1, 2, 



(125) 
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or, explicitly, 



-huj 3 < (o) |i3|p,g,s,^ >= 



rhuj 3 



4^ = ^3 < (j)^\^l\p,q,S,e >= frw 3 (-2IQM + S + 



2£QM + s 

r 



As before, we choose the JQ's in such a way that |0, 0, > is a solution, 

K t + K 2 -K 3 -^ 3 = 
K 1 -K 2 + 2K 3 + ^lu 3 = 0. 
We now require that the states have vanishing eigenvalues, 

Ku)i I p 



4 0) + Ae« 



M P A P 



+ Xe 



(i) 



iV 

M P P 



0. 



(126) 
(127) 

(128) 
(129) 

(130) 
(131) 



For these equations to have a solution with a common value of A, we need to choose carefully the states < p, q, s, £\, 
i.e, we must have q = —3Nb,p = Mb, and the solution for A is, 



A~-^ 



P£ QMl 



(132) 



so we see we completely reproduce the expansion of the exact solution we found in the previous subsection. As in the 
previous example, there a re m any eigenstates one can start from in the zeroth order theory that yield a solution for 
the same A. The formula (132) is invariant if one multiplies b, s, I times an arbitrary integer. 

We therefore see that indeed one needs extra conditions in order to have both constraints vanish on a given value 
of the perturbative parameter. But the conditions just imply that one starts with different zeroth order eigenvalues 
for the states in both constraints. The situation is schematized in the figure. 



Eigenvalues 



Eigenvalue of constraint 1 
Eigenvalue of constraint 2 




FIG. 2. A schematic diagram of how the constraints are solved. One starts from the eigenvalues of the unperturbed constraints 
of the system. The perturbation breaks the degeneracy of the eigenvalues and changes linearly the eigenvalues. When the 
corrected eigenvalues cross the zero value simultaneously for both constraints, one has a permitted value of A. 

Finally, if one sought to compute the first order correction of the quantum state, one get two equations for it, one 
per each constraint. If one studies them in detail, one finds that they admit the same solution. 
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VI. SKETCHING THE FULL QUANTUM THEORY 



We do not have a detailed discussion of the full general relativity quantum theory ready, but nevertheless one can 
sketch how the procedure would be. The calculations in this section should only be taken as a guide, and will be 
highly formal in nature. The details are very important in determining the final quantum theory, therefore we can 
only refer to these calculations as a "sketch" of what can be achieved. We will see simultaneous solutions to the 
diffeomorphism and Hamiltonian constrain. In order to achieve this we start with a basis consistent of equivalence 
classes of spin nets under subgroup of diffeomorphisms that keeps the vertices of the spin net fixed, 

\s m >=J2 D [v]\s>, (133) 

where the are the diffeomorphisms that leave the vertices v of the spin network fixed. This is a space very similar 
to the "habitat" considered by Lewandowski and Marolf . The action of both the diffeomorphism constraint and 
Thiemann's Hamiltonian is well defined in such a space. We can therefore proceed to study the problem of eigenvalues 
and eigenvectors to first order in perturbation theory, 

(< V^ (0) | + A" 1 < V^ (1) |) (V(M) + A-^(M)) = (< V> (0) | + A" 1 < |) (e (0) (M ) + A _1 e (1) (A/)) , (134) 

Considering the equation at zeroth order we get, 

< V (0) \V{M) = e(°) (M) < V> (0) | (135) 

the solution of this equation is given by the eigenvalues and eigenvectors of the determinant of the metric. Here 
we have to deal with the fact that this is a field theory. One can smear the determinant of the metric with a 
function M{x). The resulting operator is closely related to the well-understood JI|,|o| volume operator, we call it 
"smeared volume" and denote it as V(M). One will essentially get the volume associated with each vertex, which is 
a finite quantity, multiplied times the value of the smearing function at the vertex. Eigenstates will be reasonably 
easy to find. They will correspond to combinations of spin networks that yield the same total for the sum of the 
volumes of each vertex times the smearing function valued at the vertex. That quantity, in turn will be the eigenvalue 
e(°)(M) = J2v V( v i)M(vi). We can schematically write these states as < ipi^l — J2 S C(s) v < { s }v.e°(M))\ where 
we have a sum of states such that the all have the same value of V(M) for a given M. There are many states for 
each value of e^(M), to illustrate this, consider a spin network with a given number of vertices and a given value 
of e^ ) (M) and replace two lines in it by the same two lines with some knotting in between them. The value of the 
eigenvalue will not change, but it will be a different spin network. 

Let us now consider the first order corrections. As we did in the quantum mechanical examples, we will project the 
perturbative equation on the space of states < Sv(m)\i 

< ^\V{M)\s v[M) > + < ^% {M) \H(M)\s v{M) >= e«(Af) < ^l (M) \s v(M ) > +e (0) (Af) < ^% V{M) > . (136) 

The first and last term in both members cancels out, since V(M)\s V (m) >= ^°\M)\s V (m) >■ One is therefore left 
with, 

< ip% (M) \H(M)\s v(M) >= e«(M) < ^% V(M) > . (137) 

One can use this equation to explicitly compute eW(Af), for a given proposal for the Hamiltonian constraint. One can 
now proceed to make the "energy" vanish at this order of perturbation theory, therefore determining the (quantized) 
value of the cosmological constant and the possible values of the energy of the zeroth order state e°(M). 

To obtain the correction to the state, we need to project the above equation on states that are orthogonal to the 
eigenstates with eigenvalue e^(M). Let us call those states \sy(M) >> that is, states with different values of the 
volume at the vertices. We therefore assign the volume at the vertices as a measure of orthogonality. We therefore 
write, 

< 1>W\V(M)\a v , (M) > +A- 1 < i;% m \H(M)\s v , {M) >= e^(M) < ^% {m) \s v >{m) > +A-V°>(Af) < ^\s v > ( m) >, 

(138) 

Now, since V{M)\sy(M) >= e '^°\M) s v'(M) >, we can write, 
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^ ^ <^% {M) \H{M)\ S y, {M) > 

< ip v >\s v ,( M ) >= j^—^ — . (139) 

^°)(M)-c(°)(M) 



and this equation would determine the first order correction to the quantum state. This would be only strictly true if 
the state of volume V(M) were non-degenerate. The degeneracy implies that there might be non- vanishing projections 
< ip^\s V (M) > for vectors that belong to the space of states with eigenvalue e^°\M). It is well known in degenerate 
perturbation theory that higher order calculations determine these components. 

The dependence on M of the last equation might appear as surprising, since the first order quantum state should 
be M independent. Here we can draw on our experience on the system with two constraints studied in last section, 
which shows that a single correction appears no matter which Hamiltonian one uses to compute the correction. In 
the gravity case this would correspond to getting the same solution for different M's. Another way to see that the 
equation is M-indepcndent is to notice that one can characterize the action of the Hamiltonian independently at each 
vertex and in such calculations the role of M is that of a constant overall factor that drops out of the equations and 
one is left, as in the case of the oscillators, with a system of equations. 

We have therefore constructed states such that, 

<^|(A)^<V4° ) I+A- 1 <^ 1) I (140) 

such that, 

< V„(A)|#(M,A) - 0(A- 2 )VM. (141) 

One is however, interested in states that are genuinely invariant under diffeomorphisms, whereas the above states 
are only invariant under diffeomorphisms that leave the vertices of the spin network fixed. The kind of states we look 
for are, 

< V(A)| = < lM A )P- (142) 
v 

such that 

< 1>(A)\H(M, A) = 0(A~ 2 ). (143) 

We shall show that with the above defined states this is indeed the case. This follows from the Poisson algebra of 
diffcomorphism and Hamiltonian constraints, 

VH(M, A)X> _1 = H(VM, A). (144) 

Therefore, 

< i[)(A)\H(M, A) = < ipv(A)\'DH(M, A) = ^ < il> v (A)\H(DM, A)T> = 0(A~ 2 ), (145) 
v v 

which is what we wished to show. 

The details of these calculations can only be filled in with a given prescription for a Hamiltonian constraint and a 
definite space of states to operate upon. This will require much more detailed work than is appropriate for this paper. 
One can hint at what would happen, for instance, if one uses Thiemann's Hamiltonian. Let us start with trivalent 
spin networks. In such a context V(M) vanishes identically and the perturbative method has nothing to add: one 
is left without higher order corrections and the solution to the problem is reduced to determining which states are 
annihilated by the Hamiltonian constraints. To get something non-trivial, one needs four-valent vertices. Calculations 
on four valent vertices are readily possible, but cumbersome. In this context one could sharpen the determination 
of eigenstates, etc. Very schematically, one would see that the correction to the states that one gets depend on the 
action of the Hamiltonian constraint at each vertex. The states appear different from the ones obtained via "group 
averaging" that involve acting repeatedly with the Hamiltonian. However, since the action of the Hamiltonian is 
"local" (i.e. it does not add connections among vertices, but "dresses up" each vertex), the kind of states one gets 
have similar locality properties to those obtained with the "group averaging" procedure. In fact, if one works at higher 
order, at each order one acts with a Hamiltonian, so one can see that by going to higher orders, one is recovering 
something similar to the "group averaging" . 

An interesting point in trying to draw an analogy between quantum gravity and the simple models discussed in this 
paper is if the perturbation Hamiltonian has a discrete spectrum or not. We know that the zeroth order Hamiltonian 
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(the smeared volume) has a discrete spectrum. At the moment, it is knot known if either the Thiemann or the Vassiliev 
Hamiltonians have a discrete spectrum. However, we have shown in this paper that the perturbative method works 
well in the case in which the perturbative Hamiltonian has a continuous spectrum (the case of the coupled oscillators). 
It is well known that the perturbative approach works well when the perturbing Hamiltonian has a discrete spectrum 
(for instance if one considers an atom in a magnetic field with constant energy). 

It is clear that this is just the beginning of a discussion of the four dimensional quantum gravity case. It is interesting 
to notice, however, how one can quickly gain intuition as to the kind of states one would get and how they depend 
on the action of the Hamiltonian without doing the explicit computations. Since the latter must involve four valent 
vertices and therefore are extensive in nature, having a quick way to intuitively try out proposals for the Hamiltonian 
is a great asset. 



VII. CONCLUSIONS 



It is obvious that the complexities that one expects in quantum gravity (in particular the infinite dimensional 
nature of the problem) cannot be really mimicked by finite dimensional systems. We believe, however, that the 
finite dimensional models we analyze in this paper help dissipate some of the most elementary skepticisms about the 
possibilities of the approach we are proposing. Namely: 1) that the method seems to simplistic to deal with chaotic, 
possibly pathological Hamiltonian systems as one expects full general relativity to contain in certain regimes; 2) that 
the method is only circumscribed to a range of values of the cosmological constant of no physical relevance; 3) that 
the application in the quantum domain of the technique is possibly inconsistent. 

Applied to the full quantum theory, our method appears to reduce the problem of finding quantum states to a set 
of well defined, albeit complicated spin network calculations. Moreover, it appears to yield a quick intuitive handle 
on the properties of the quantum states of possible proposals for Hamiltonian constraints. We are currently exploring 
further the application of the approximation technique in full quantum gravity, as sketched in the last section. 
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